Fermi velocity renormalization and the excitonic insulator in graphene 
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We compute the Fermi velocity of tiie Dirac quasiparticles in clean graphene at the charge neutrality point 
for strong Coulomb coupling a^. We perform a Lattice Monte Carlo calculation within the low -energy Dirac 
theory, which includes an instantaneous, long-range Coulomb interaction. This method is non-perturbative and 
takes full account of quantum fluctuations. We find a finite renormalized Fermi velocity Vpn > Vp, where 
Vp ~ c/300. We introduce the critical Fermi velocity renormalization v^ = Vpft{ctgc)/'>^F' where a^^ is the 
critical coupling for the semimetal-insulator transition due to excitonic pair formation. We compare our results 
with empirical studies of interaction-induced spectral changes in graphene. We find v^. ~ 3.3, which should 
be contrasted with Vp^/vp ~ 2 — 3 for ultra-clean suspended graphene and Vp^/vp ~ 1.2 for graphene on 
a boron nitride substrate. Our results are consistent with the non-observation of insulating states in suspended 
graphene in the absence of an external magnetic field. We also discuss the dynamical critical exponent z. 
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Graphene, a two-dimensional membrane of carbon atoms 
with unique electronic properties, is characterized by a low- 
energy spectrum of Dirac quasiparticles, with a Fermi veloc- 
ity Vp ~ 1/300 of the speed of light in vacuum ||jj,|2|]. As 
the strength of the Coulomb interaction between the quasi- 
particles is controlled by a = e^ /{Airevp) ~ 2.2/e, the 
role of interactions can be enhanced to the point that graphene 
may resemble Quantum Electrodynamics (QED) in a strongly 
coupled regime HSfl. In particular, the unscreened, long-range 
Coulomb interaction in graphene leads to non-trivial veloc- 
ity renormalization effects. These have been extensively ana- 
lyzed in the weak-coupling regime Wl, and at strong coupling 
in the random-phase approximation ||5l|6|]. At weak coupling, 
a logarithmic running iJ^(n)/wp(ng) ~ 1 + (a„/4) ln(nQ/n) 
with carrier density n is found, thus the renormalized Fermi 
velocity Vppj is expected to become large in the vicinity of 
the Dirac point (n = 0). On the experimental side, loga- 
rithmic velocity renormalization has been detected in ultra- 
clean suspended graphene ||7|], on boron nitride (BN) sub- 
strates |8], and in ARPES measurements of quasi-freestanding 
graphene on silicon carbide (SiC) H. These experiments find 
Vpj^/vp > 1 at the Dirac point, with the largest renormaliza- 
tion Vppi/vp ~ 2 — 3 in suspended graphene, where e = 1. 

In addition to the renormalization of Vp, strong electron- 
electron interactions may trigger a semimetal-insulator tran- 
sition due to the excitonic pairing of quasiparticles and holes 
above a critical coupling a„^. Notably, spontaneous gap gen- 
eration in graphene has been studied using the Schwinger- 
Dyson (or gap) equations llioll . and within the Lattice Monte 
Carlo (LMC) framework, which is non-perturbative and takes 
full account of quantum fluctuations. In particular, LMC is 
noted for its central role in the study of dynamical mass gen- 
eration (or spontaneous chiral symmetry breaking) in gauge 
field theories, such as QED and Quantum Chromodynamics 
(QCD), see Ref. flllll and references therein. In the case of 
graphene, LMC has been applied to the Dirac theory using a 
contact interaction related to the Thirring model [:12], a long- 



range Coulomb interaction ||I3LII4|1 . and to the tight-binding 
description using a Hubbard interaction ||l5||. The LMC cal- 
culation of Ref. 11311 found a ~ 1.11 ± 0.06, which sug- 
gests that the insulating phase might be observable in ultra- 
clean suspended graphene samples, as a < 2.2. So far, the 
semimetal-insulator transition has not been detected, although 
the observed Vpp^/vp ~ 2 — 3 in suspended graphene is in- 
dicative of interaction-induced spectral changes 131 17|]. 

Our objective is to compute Vpj^{ag)/vp at the Dirac point, 
and to contrast this with the electronic phase diagram. Specifi- 
cally, we shall consider the "critical velocity renormalization" 
{a )/vp, i.e. the velocity renormalization at the 



critical coupling a . Then, the semimetal-insulator transition 
should only be observable when the velocity renormalization 
exceeds v^ for n = 0. A clear advantage of v^ over a is that 
the former is a renormalized (physical) quantity which can be 
directly compared with a wide range of experiments. 

The LMC treatment of graphene uses the linearized low- 
energy Hamiltonian jm llTIl with an instantaneous Coulomb 
interaction. This gives the Euclidean (continuum) action 
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where g^ = e^/e with e = (1 + k)/2 for a substrate with 
dielectric constant k, ifj^ is a four-component Dirac field in 2+ 
1 dimensions with ip = i/'^To' ^o ^^ ^^^ gauge (photon) field 
in 3 + 1 dimensions, and N s = 2 for a graphene monolayer. 
The Dirac operator is 



D[Aa] = 7o(ao + iAo) +vfJ2 ^kdk 



'Oi 



(2) 
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where the Dirac matrices 7 satisfy the Euclidean Clifford al- 
gebra {7m, 7i/} = 2(5„j,, and m^ is the bare fermion mass, 
which provides an infrared regulator for modes that would 
otherwise be massless when the \5{2N A chiral symmetry of 
Eq. (dJ is spontaneously broken. 



The lattice version of Eq. ([T]l is preferably formulated 
in terms of "staggered" fermions llisll , i.e. one-component 
Grassmann variables x, x, which partially retain the U(2A^,) 
chiral symmetry of Eq. ([T) at finite lattice spacing. The full 
symmetry is recovered in the continuum limit. Every other 
lattice site in the staggered formulation is identified with a 
space-time degree of freedom, and the resulting "hypercubes" 
of lattice sites are treated as internal degrees of freedom. The 
four-component Dirac spinors are reconstructed in the contin- 
uum limit as appropriately chosen combinations of the stag- 
gered fermions in each hypercube. Staggered fermions have 
been successfully applied to theories ranging from QED3 lll9ll 
and QED4 liolHil to the Thirring model EH. In QCD, 
staggered fermions have been developed into highly improved 
lattice actions, which systematically remove the discretization 
errors, see Ref. 11 ill for an overview. The application of stag- 
gered fermions to Eq. ([T]i is detailed in Ref. [11311. and a peda- 
gogical introduction can be found in Ref. 112511 . 

The fermionic part of Eq. dU is given for Nr = 2 in terms 
of staggered fermions as J2n,m Xn -f'^n.mi^o] Xm- where n = 
(riQ, 'ni,n2) = {t, x, y) and m denote lattice sites on a 2 + 1 
dimensional fermion brane. The staggered Dirac operator is 
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where ?y^ = (—1)"", rj^ = (— i)"o+"i^ with e^ a unit vector 
in lattice direction i. The invariance of Eq. ([TJ under spatially 
uniform, time-dependent gauge transformations is retained by 
coupling the fermions to the lattice gauge field 9q via "gauge 
links" Uq = cxp(z6'q). Our LMC calculations are performed 
for 1)^ = 1 (thus g^ -^ 9^/vp) and we take a/a^ ~ 1, where 
a = Cj is the temporal lattice spacing, and therefore A = 1. 
However, at non-zero q;„ we have A^ = Vpj^ {a/a.^)pi, from 
which the velocity renormalization Vp^/vp can be obtained 
once the lattice spacing asymmetry {a/a^)p, is known. We 
emphasize that our cubic space-time lattice should not be un- 
derstood in terms of a square tight-binding theory. The low- 
energy physics of graphene described by Eq. ([1} is recovered 
in the continuum limit a — ?> 0, taken along lines of constant 
physics (as a function of the bare parameters) in the direction 
of increasing correlation length lll3[ |21 1 . 

We compute the renormalized A^j and rrip, from the stag- 
gered fermion propagator 

Cf{t,x,y) = {x{t,x,y)x{t„,Xo,yo)) = {K^^nJ, (4) 

on an N^ x N^ space-time lattice with N^ j/4 integer Here 
Hq is an arbitrary point of reference, and the brackets denote 
an average over "ensembles" of gauge field configurations, 
obtained as a function of /3 = Vp/g'^ = l/(47ra ) and ttiq 
with the Hybrid Monte Carlo (HMC) algorithm U^- Once 
CJt, X, y) is computed, we form the temporal correlator 

Cft{t,Pi,P2) =^exp{-ip-x)Cj:{t,x,y), (5) 
x,y 
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FIG. I: Calculation of A^ as a function of /3 = l/(47rag) and mg 
from the correlators Cj^ and Cf^ on different lattices. Statistical 
errors are comparable to the size of the symbols. 



for momenta p^ = 2im/N^, with n = —N^/A, . . . , N^/4, 
timeslices t ~ 0, . . . ,Nf~l, and summation over even lattice 
sites, defined by (— l)*+^+y = 1. The temporal correlator CS 
with renormalized tti^ and A^ has been given in Ref. 12111 . 
The expression for C5 with wave function renormalization 
Zjj. is the "sawtooth" function (for a = 1) 



C/t(t,Pl,P2) = ^i?G't(pi,P2), 



(6) 



fort = 0,2,...,A^t 

Cft{t,Pl,P2) 
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cxp(iBo)Gt+i(pi,P2) 



exp(-iBo)G't_i(Pi,P2) 



(7) 



for t = 1,3, . . . ,N^ — 1, with anti-periodic boundary condi- 
tions. The function Gf{pi,P2) is given by 

GM,P2) = N[cos{B,Nj2)cosh{ti,Nj2) 

X cosiB^t*) sinh(/Xti*) + sm{BQNj2) smh{fj.tNj2) 
X sm{B„t*) coship^t*) + i cos{B()Nj2) cosh{fitNj2) 
X sm{B„t*) smh{fi^t*) ~i&\n{BQNj2) fim\v{^i^N J 2) 
X cosiB^t*) cosh{fj,^t*)] , (8) 



where t* = Nj2-t, 
N ^ . ^7f [cosh^(A.,iV,/2) - sin^B,Nj2)y' , (9) 

and the dispersion relation 

sinh^(Mt) = TO^ + A^ siii2(pi) + A^ sin^ip^)- (10) 

The expression for C?^ includes a constant "background field" 
Bq = {9q) in CS, as (0q) 7^ in a finite volume. The back- 
ground field Bq is roughly bounded by ±n/Nf 12111 . Note that 
the imaginary part of CS vanishes in the limit Bq — > 0. 

We also consider spatial correlations, to evaluate the effects 
of anisotropics. We define the spatial fermion correlator 



Cf^{x,ui,P2) = ^ exp{-ip ■ x) Cf{t, X, y), 



(11) 



for uj = 27r(ri — l/2)/Nf (due to the anti-periodic tempo- 
ral boundary conditions), n = —Nf/A, . . . , N^/A, spaceslices 
X = 0, . . . , N^ — 1, and summation over even lattice sites. 
The expression for C?^ can be inferred from C5. The func- 
tion G^(a;,p2) is obtained from Eq. dHJ by first exchanging 
sin o cos, followed by the substitutions t ^ x, Nf —^ N^, 
Hi —!' /ij. and Bq —^ 0. These also apply to Eq. O and (|9]l, in 
addition to m^ — > rrij^/Xj^ in Eq. dTji and m^^ — ?■ irij^/Xj^ in 
Eq. (|9]l. The dispersion relation is 



sinh^ifi.^ 



j/A' 



jj + sm'{u; + B„)/Xjf + sm'ip^). (12) 



Unlike Eqs. ^ and (|7]i, C?^ is real-valued and satisfies peri- 
odic boundary conditions. 

In Fig. [T] we show A^ as obtained from a chi-square fit 
of Eqs. © and ^ to LMC data. While Eq. ([B is gauge in- 
variant, the fermion correlators Cf^. and C^^ are not, and thus 
a gauge fixing condition is imposed on each configuration in 
order to obtain stable results. For C^j, a correlated fit is per- 
formed for all t, pi , P2, and in the case of C,^ for all x,uj,P2- 
The fitted parameters are m^, A^, Bq and Zj^. Our lattices 
have Nf. = N^, and the length of the "bulk" dimension (where 
only the photons propagate) is denoted N^^. We use the nota- 
tion TVj X Nf^, and simply N^ when N^ = N^ = N^. On the 
28^ X 8 lattice, LMC data is available for (inverse) lattice cou- 
plings 5.0 < f3~^ < 15.0, and for bare quasiparticle masses 
0.0025 < ttiq < 0.020, with slightly more restrictive data sets 
on the 32'^ x 12 and 32^ lattices. We find that A^ increases 
monotonically as a function of a„ from the non-interacting 
value of unity, with no appreciable differences between A^ as 
obtained from Cff and Cf^. We find the dependence on ttiq 
to be almost negligible. Finite size effects for A^^ are small, 
and the fitted values of Bq agree closely with {9q). 

In Fig. 121 we show the physical quasiparticle mass nij^ as a 
function of /3 and ttiq, with emphasis on asymmetries between 
the temporal and spatial correlations, and on finite size effects. 
As Eq. ([T]i treats space and time asymmetrically, the spatial 
and temporal correlation lengths ^ may exhibit unequal criti- 
cal scaling, such that ^j, oc |/3 — /3^|~'^= and ^^ oc |/3 — f3^\~^t. 
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FIG. 2: Calculation of the physical quasiparticle mass niji as a func- 
tion of j3^^ = A-KOg and mo (for a = 1). Left panel: Calculation of 
TTifl from temporal (colored connected symbols) and spatial (black 
unconnected symbols) correlations on a 28'' x 8 lattice. Right panel: 
Calculation of rriji from temporal correlations on a 28'^ x 8 lattice 
(colored connected symbols) and on a 32^ x 12 lattice (black uncon- 
nected symbols). All lines are intended as a guide to the eye, and 
statistical errors are comparable to the size of the symbols. 



The dynamical critical exponent z = Vi/Vg is an important 
characteristic of a quantum critical point (QCP), and implies 
that the dispersion relation is modified to E oc p^. At large 
Nr, Ref. 16(1 predicted z ~ 0.8 for graphene in the strong- 
coupling limit. However, arguments have also been put for- 
ward which indicate z = 1 for a QCP with d < 4 in theories 
with a long-range Coulomb interaction ll27ll . From Fig. |2] we 
find that the values of rrijj obtained from Cf^ and Cf^ agree 
very closely for /3^^ < 11.0, which is consistent with z ~ 1. 
We also find no sign of non-linear dispersion. A more accu- 
rate analysis is possible following Ref. II12II . in terms of the 
equation of state (EOS) 



Aif3 



/3j-r" 



-1)/"^* 
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(13) 



for m^ computed from C^j, where 5 and /?„, are critical ex- 
ponents characterizing the QCP at /? = /?^. An analogous 
EOS with f( — )• Vg can be given for mp, as obtained from 




FIG. 3: Plot of vpji{ag)/vp for asymmetry {a/a^)^ ~ 1. The 
critical velocity renormalization is v,. ~ 3.3 at a^^ ~ 1.1, and the 
predicted insulating phase occurs for Qg > Ogc (gray shaded area). 
We find a linear dependence of Vp^/vp on a^ up to a^ ~ 0.5 (solid 
black line). Note that vpii{0)/vp = 1. Horizontal bands indicate 
empirical data on vpji/vp from Ref. [3] (suspended graphene) and 
Ref. [3] (graphene on a BN substrate). 



corresponds to a ~ 0.7 — 0.9, although it should be kept 
in mind that a„ is an unrenormalized parameter. While the 
electron charge is believed ll28ll not to be renormalized, the 



calculation of e^ in graphene remains a challenge for LMC. 

The recent study of the velocity renormalization in extrinsic 
graphene (?7 ^ 0) by Ref. 12911 has predicted a dramatic break- 
down of the quasiparticle picture, which should be observable 
in graphene on an Si02 substrate (a ~ 0.8), and dominate 
the spectral properties in suspended graphene {a — 2.2). 
While we find no breakdown of the description in terms of a 
renormalized quasiparticle propagator, our LMC data extends 
only to a ~ 1.2, thus we cannot at present determine whether 
the quasiparticle picture remains valid at larger couplings. 
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for instructive discussions, and Jan Bsaisou, Dean Lee and 
Ulf-G. MeiBner for a careful reading of the manuscript. T. L. 
acknowledges financial support from the Helmholtz Associa- 
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by an allocation of computing time from the Ohio Supercom- 
puter Center (OSC). 



Cf^. We also find from Fig.|2]that finite size effects are under 
control for f3^^ < 9.0, but not for smaller /3 (stronger cou- 
pling), especially in the region of the phase diagram where a 
dynamically generated gap exists, i.e. liin„i .^0 "^_r t^ *-*. In 
principle, Eq. (fTsT i can be used to compute z and rrij^ in the 
limit niQ -^ 0, i.e. the gap in the quasiparticle spectrum. The 
lattice spacing asymmetry (0/0^.)^. could then be inferred by 
measuring the gap in terms of temporal and spatial correla- 
tions. For reliable conclusions, such an EOS analysis should 
be combined with an extrapolation to infinite volume, similar 
to that of Ref. 12211 for QED4. For the present analysis, we 
note from Fig. |2]that the values of irij^ computed from Cf^ 
and Cf^ differ by no more than ~ 10% at the smallest values 
of /?. Assuming this trend persists in the limits of infinite vol- 
ume and vanishing TOq, we find 1.0 < (a/a^)jj. < l.loverthe 
range of a studied. In the absence of substantial indications 
for {a/a^)p ^ 1, we take A^^ as a measure of Vp^Jvp. 

In Fig. [3] we summarize our LMC results for Vpj^/vp as 
a function of a = 1/{4tt/3), and compare with available 
experimental data. We find that Vpj^ increases linearly with 
a from the non-interacting value Vp up to a ~ 0.5, above 
which the increase becomes more rapid. At the predicted crit- 
ical coupling Ug^ ~ 1.1, we find v^ = Vpn{agJ/vp ~ 3.3 
for the critical velocity renormalization. Since all of the em- 
pirical Vpj^/vp fall short of v^, we find a plausible explana- 
tion for the non-observation of exci tonic insulating phases in 
graphene monolayers. Nevertheless, the result of Ref. [7]] is 
tantalizingly close to v^, which suggests that further refine- 
ments in the quality of suspended graphene may suffice to 
trigger the excitonic instability. From Fig. |3] we find that 
2 < Vpn/vp < 3 reported by Ref. [7]] in suspended graphene 
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